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PACS 04 . 62 . +v - Quantum fields in curved spacetime 

Abstract. - We find out classical particles, starting from Dirac quantum fields on a curved space- 
time, by an eikonal approximation and a localization hypothesis for amplitudes. We recover the 
results by Mathisson-Papapetrou, hence establishing a fundamental correspondence between the 
coupling of classical and quantum spinning particles with the gravitational field. 



I- Introduction. — Mathisson [1] and Papapetrou [2] 
described the behaviour of a spinning particle in general 
relativity (GR). The work of Papapetrou is particularly 
worth noting. He performed a multipole expansion around 
the worldlinc. The motion turns out to be non-geodesic, 
because of the interaction between spin and curvature, 
while the evolution of the spin consists in a precession 
around the generalized momentum. 

Since the Mathisson-Papapetrou system does not repro- 
duce the classical force for the spin-orbit interaction [3], 
Pomeransky and Khriplovich [4] adopted a different pro- 
cedure. Given a distribution of energy momentum T^ v in 
a background metric g^ v , they imposed the conservation 
of the spin tensor along the trajectory and a coupling term 
of the form WT^, 

— Qy,v j Vfj.^ being Minkowsky 

metric. 

The dynamical equations obtained within these two frame- 
works are different. 

On the quantum level, Audrescht [5] constructed the 
semi-classical limit of the Dirac equation on a curved 
space-time using the eikonal approximation (the case with 
a non- vanishing torsion is presented in [6]). If the classical 
4-momentum was identified with the convection current, 
the equations dexribing the trajectory and spin evolution 
turned out to be the Papapetrou equations. However, the 
transition to the classical picture in terms of a localised 
worldlinc was not discussed. 

Recently, Silenko and Teryaev [7] found, applying a Foldy- 
Wouthuysen transformation [8] on the Hamiltonian, that 



the equations describing the evolution of the spin and of 
the momentum associated to a semi-classical spinor in a 
weak gravitational field coincide with the Pomeransky- 
Khriplovich equations. 

In this work, the classical-quantum spinning particles 
comparison is made only when a proper classical descrip- 
tion is inferred from quantum theory. This implies that 
the quantum field must be localized, such that a trajectory 
can be properly defined. In this respect, we face a second 
order approach to the semi-classical limit of Dirac parti- 
cles moving on a curved space-time by an eikonal approxi- 
mation and a localization hypothesis along the world-line. 
The latter is obtained to the lowest order of a multi-pole 
expansion la Papapetrou (single-pole), such that by spa- 
tial integration one deals with quantities having support 
on the world-line. Within this scheme we are able to de- 
rive a dispersion relation from the squared Dirac equation 
such that general covariance is manifest. No restriction to 
weak gravitational fields is imposed. Finally we retain the 
leading order in the semi-classical expansion, which means 
the first order in the spin, and we get a set of equations 
which coincide with the Mathisson-Papapetrou ones. 
For a different approach to recover Papapetrou-like equa- 
tions from the semi-classical dynamics of Dirac particles, 
in the context of using Grassmanian variables, see [9], 

II- Papapetrou formulation. In GR a moving 
body is described by an energy-momentum tensor having 
support on the world-tube r. In order to determine the 
dynamics of the body and of the gravitational field, one 
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should solve the corresponding set of Einstein equations. 
However, the solution of this system of partial differential 
equations is a rather difficult task, therefore a very useful 
approximation consists in neglecting the back-reaction of 
the body, by assigning a fixed space-time geometry. 
This way the full dynamical information is encoded in the 
energy-momentum conservation V M T M,y = 0. Inside r one 
identifies a curve X fJ -(s) and calculates momenta of the 
body [2]. 

The multi-pole expansion consists in an expansion around 
X^^s). At the zero-order one retains only the first mo- 
ment, which means that the gradients of the gravitational 
field arc neglected over the whole spatial extension of the 
body. This way, X fJl (s) results in a geodesic curve, as 
expected. At the next order of approximation (the pole- 
dipole case) gradients of the Christoffel connections are 
no longer neglected. Hence the following quantity (spin 
tensor) 

= J 6x»T v0 dT - J 5x v T^°dT (1) 

is non-vanishing and it enters into equations of motion 
giving 

JL PP — I ft PQP&TT" 
Ds 2 pav u u 

D-S» v = P^U V - P V U^ , (2) 

pi = mU V. _ DS^ Uv 

m and being the mass and the 4-velocity of the body, 
respectively, while R pa ^ is the Riemann tensor. 
The main conclusions arising in this approach can be sum- 
marized as 

• the dynamics is described by the generalized momen- 
tum P^, which in general is not aligned with the 4- 
velocity U^; 

• the spin tensor motion consists of a precession around 
the 4-velocity; 

• the trajectory is not a geodesic, since there is an in- 
teraction term between the gravitational field and the 
spin tensor (right side of the first equation in (2)). 

To provide a solution of the system (2) three conditions 
must be added, because there are thirteen unknown quan- 
tities and ten equations only. In the literature one finds 
three possible choices, 

S H0 = Qj sl *u Uv = Q ; S nvp v = Qj 

the Corinaldesi-Papapetrou [10], the Pirani [11] and the 
Tulczyjew [12, 13] consistency conditions, respectively. 

Ill- On the localization hypothesis. The de- 
scription of a fundamental particle in terms of a classical 
trajectory around a picked packet is affected by the spread 
of the wave function, in view of its quantum nature. 
We are going to outline that even though this process 
occurs, nevertheless a one-particle description can be in- 
ferred on sufficiently large macroscopic scales for certain 



initial conditions. 

The key point of our treatment consists in the existence 
of two length scales. The first one gives the initial wave 
packet spread and we denote this scale with A = a\ c , 
A c being the Compton scale, while a > 1. The second 
length scale we deal with is the radius of curvature of the 
space-time manifold, i.e. the one associated with the grav- 
itational field. We indicate it by L. In view of applying 
the theory of quantum fields on a curved space-time, we 
require L » a\ c . This condition is necessary to recover 
a proper notion of particle and when violated we enter into 
the quantum gravity regime. 

Given a time-like Killing vector, the energy of the par- 
ticle can be defined by a dispersion relation of the form 

to = \f w + g ij kikj. 

Starting with a wave-packet centered around the wave vec- 
tor k — (fco,0,0), the corresponding wave function reads 
as follows 

<f,(x,t) oc j d 3 ke-* 2 (kl ~ k < +k22+k K ^-^ W*, (4) 

1/ A being the spread in the momentum space (we take 
it as isotropic). As far as a Taylor expansion around ko is 
concerned, by neglecting terms of the (ki — (fco)i) 3 order, 
the wave function turns out to be localized around x = 
(vct,0,0), i.e. [14] 

-, (a: 1 -vet) 2 (x 2 ) 2 + ( a: 3 ) 2 

<t>{x,t) = tjz — e , (5) 

v = ko/oj. Therefore, the expectation value describes the 
trajectory of a classical free-particle. However, this is not 
enough to recover a classical particle dynamics, since the 
wave packet spreads, thus becoming less and less localized 
in space. This spread is due to the monotonic increase 
of deviations a p and a t along the directions parallel and 
orthogonal to the trajectory, respectively. 
Their expressions can be evaluated, finding 

/ i i? 

<tyt(*) = \/ A " + X P/tjIy c2t2 i X P = ^2 A * = L 

(6) 

As a consequence of this spread, soon or later the classical 
description is no longer available and a trajectory cannot 
be inferred. The point at which this failure happens de- 
pends on the minimum distance the detector can probe 
and hence can be made arbitrarily large by using a less 
accurate device. 

However, our aim is to make a comparison between the 
behavior of classical and quantum-like spinning objects, 
as far as their interaction with the gravitational field is 
concerned. Hence we are interested in a formal definition 
of a trajectory, which depends only on properties of the 
beam. We will outline in the Discussion the implications 
of this choice on the experimental level. 
For all the speculations above, we take the distance trav- 
eled by the expectation value as the characteristic length 
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scale with which the spread has to be compared. There- 
fore, we set the following condition to yield a one particle 
dynamics: 

vet » a(t) - A. (7) 

The behavior of x 1 = vet (solid line) and a(t) — X are given 
by curves in figure 1. 




Figure 1: The solid line gives the expectation value of the posi- 
tion, while the dotted, dash-dotted and dashed ones correspond 
to the variation of the spread, i.e. a(t) — A, in cases av < 1, 
av = 1 and av > 1, respectively. It is worth noting how for 
av > 1 the spread always keeps smaller than the proper time 
along the world-line of the wave packet center. 

It is worth noting that while for av < 1 a one-particle 
description drops down after a distance vt* w (aw) 2 A < 
A c , for av > 1 the spread is always much smaller than the 
distance. Hence, the localization holds until the quantum 
gravity regime is approached. 

Therefore, in view of speculations above we deal with wave 
packets for which 



a > 1/v = 1 + 



1/2 



(8) 



This condition stands for massless particles, while the ini- 
tial spread must be larger as the mass increases. How- 
ever, in view of the great difference between the Compton 
scale and that one characteristic of macroscopic phenom- 
ena, there is a large set of velocities available, without 
affecting the notion of particle at the initial time. 
We emphasize that this sort of localization allows us to 
treat particles as being point-like and to neglect spatial 
gradients in a proper reference frame (see next section) 
up to the notion of a trajectory is lost. In this approx- 
imation scheme the time variation of the spread can be 
neglected too and in what follows we will treat it as con- 
stant. Indeed this description drops down for observers 
which are able to test scales much smaller that L. 

IV- The semi-classical limit of the Dirac equa- 
tion. — The Dirac equation in a curved space-time reads 
as follows (we work in units h = c = 1 with signature 



(+,---)) 



i-fD„ V = MV, D ^ = ^ - rf^ab 1>> (9) 



where we have rewritten spinor connections in terms of 
the gravity spin connections = e V^e" and Lorentz 
group generators S a b = \ [j a > lb] ■ If we square the expres- 
sion (9), by using relation \D ll ,D v \ij) 
obtain 



we 



(10) 



R being the scalar curvature, i.e. R = Rp V = 
g»» ■ gP° R jlpva and the Ricci tensor Rp V and the Riemann 
tensor Rp Vpa retain the expressions as in [17]. In what fol- 
lows, we will address the last equation as the starting point 
of our dynamical analysis. Although by passing from first- 
to second-order equations the set of solutions is enlarged, 
nevertheless in this semi-classical context this procedure 
has to be regarded as a correct one. In fact, it allows us to 
treat spinor particles as scalar ones with an H correction 
due to the spinor structure. 

Nevertheless, it is worth noting that differences exist be- 
tween our method and the one carried on in [7] , where the 
semi-classical limit is performed by a Foldy-Wouthuyscn 
transformation. 

The semi-classical limit of a quantum field dynamics 
consists of two steps: the eikonal approximation, by which 
a classical divergent phase is introduced, and the local- 
ization hypothesis for amplitudes, from which a classical 
trajectory is inferred. 

The eikonal approximation for spinors reads as follows [16] 



V> = e iS ' h u, 



(11) 



u being a spinor, for which we assume the proper local- 
ization along the world-line = X^(s), i.e. we take the 
product form 



u(s,x) = n^ =0 



1 



2ir<r, 



-e u . (12) 



We take A M (s) as the integral curve of the wave- vector 
Let us now focus on the expression 



/ 



d 3 xy/^g 



(0), 



(13) 



where the Lorentz frame co-moving with has been cho- 
sen, i.e. Up = ep°\ Up being the 4- velocity. 
We point out that the expression (13) vanishes as a con- 
sequence of squared Dirac equation (10). 
Hence, by inserting the form (11) for ip, u given by (12), 
we find 



(KpK» - K^Sp){\ + 0(X 2 )) + /i 2 = 



(14) 
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where the 0{\ 2 ) corrections come from the expansion of 
the gaussian around the world-line, while Sp reads as 



2i 



(15) 



Hence the dynamics of is obtained by acting on the 
relation (14) with the derivative operator V '„ and we have 



2 ^pafiiy 

P u = K v - S, 



o{\ 2 ) = o 



(16) 



where S^ v is the charge associated with the spin density 
of the Dirac field, whose expression reads 



2 / d 3 xVhwy u 



2u j°u 



+ 0(\ 2 ). 
(17) 

Furthermore, we outline that S^ u couples with the Rie- 
mann tensor as the classical spin tensor. Such an identifi- 
cation is enforced by the following relation 



S V »U V = 0, 



(18) 



which coincides with Pirani condition (3). 

The hypothesis of dealing with a multi-pole expansion 
is now translated into a proper condition on uq. We saw in 
the mono-pole case that the gravitational field can be ap- 
proximated as being constant over the whole spatial exten- 
sion of the body. In the same way, here uq can be treated 
as a parallel transported spinor on the space orthogonal 
to Ufj,. In this respect, the most general expression one 
can take for uq is the following one 



(19) 



v r being an arbitrary spinor, while r is the spinorial in- 
dex. The condition above is well-grounded by virtue of the 
analysis of the previous section. We do not expect uo to 
appreciate any sensible dependence on spatial coordinates, 
since in the adopted approximation scheme the full wave 
function u has a negligible amplitude outside = X^(s). 
On a physical point of view, the particle does not spread 
to much to "see" spatial gradients, but it travels so far to 
interact with the curvature along the world-line. 
This way, equations giving the dynamics can be rewritten 
as follows 



r/^v P — £ 7? 

IfPVpS^ = 
P V = K V - ns v 



U»SP a = 



(20) 



It is worth noting that from the relation (19) Pp turns 
out to be aligned with Up. Furthermore, being Sp pro- 
portional to Up, the Tulczyjew condition stands as a con- 
sequence of the relation (18), hence the trajectory is well- 
defined. 

In order to determine the dynamics we rewrite Pp = 



(n+n)Up and we obtain by the first equation of the system 
(20) 

D „ D 



a — U u H (nU u ) 

^Ds M Ds y ^' 



0, 



(21) 



which gives -jy^n = 0. Hence by a re-definition of the 
mass, the equation giving the trajectory turns out to be 
as follows 



U p V p S^ = 



(22) 



Hence the obtained dynamics is consistent with Papa- 
petrou equations (2) at the lowest order in the spin. 

Therefore, Dirac particles follow the trajectory of clas- 
sical spinning ones, whose spin tensor is given by 
(17). An interaction with the gravitational field is pre- 
dicted. We emphasize that our results differ with respect 
to [18], where equations of motion are inferred from a semi- 
classical Lagrangian density. 

The whole approximation scheme is consistent, since 
spinors can "see" only those components of the Ricmann 
tensor with at least one time-like index, i.e. those ones 
which do not contribute to the tidal forces on the spa- 
tial hypersurfaces. Such a picture has been suggested by 
results of the previous section. We have seen that the one- 
particle approximation stands as soon as at) > 1, while 
quantum gravity enters into the physical description for 
t ~ avL. Therefore, we expect that a can fixed such that 
the particle has enough time to interact with the geometry 
along its trajectory, while the space on which the packet 
spreads remains effectively flat. 

For instance, the analogous conclusion can be in- 
ferred from the classical analysis of the trajectory in a 
Schwarzschild space-time. In this respect, in [21] it is out- 
lined that the deviation from a geodesic motion is due to 
the interaction of the spin with the gravito-magnetic part 
of the metric, while tidal forces do not affect the dynamics. 

Finally, we want to stress that the assumptions made on 
the wave function (11) (i.e. the eikonal approximation and 
the localization hypothesis) are crucial to infer the above 
result. Hence the application of the Ehrcnfcst theorem 
would not be enough to obtain a Mathisson-Papapetrou 
trajectory. 

V- Discussion. — We outlined in section III that our 
results allow to make a formal comparison between the tra- 
jectories of quantum and classical spinning objects. Even 
though in principle experiments can be realized to probe 
such a model, the magnitude of the predicted effects is too 
weak to be detected. 

The main consequence of the picture outlined is that the 
particle trajectory deviates from geodesic motion. This 
deviation into a Schwarzschild space-time is discussed in 
[22] (see formula (35)), where the Pomeransky-Kriplovich 
and the Mathisson-Papapetrou results are compared. The 
main conclusion of that analysis is that different deviations 
are predicted and only in the first case the classical spin- 
orbit interaction is reproduced. 
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However, for Earth-based experiments the tidal forces are 
so weak that the associated geodesic deviation takes place 
on length scales remaining within the dimension of the 
packet. As far as particles coming from extra-galactic 
sources are concerned, even though we are supposed to 
deal with beams having a proper energy (no pair produc- 
tion), the semiclassical description is not more available 
for observers on the Earth. Another possible experimen- 
tal confirmation could come from COW-like [19,20] exper- 
iments. Since the spin couples to the curvature, trajecto- 
ries at different height are characterized by a different time 
of flight, which add a further contribution to the quantum 
interference. However, such an additional term is well be- 
low experimental uncertainty [23]. 

A different way of detecting the interaction between the 
spin and the gravitation field is the appearance of an effec- 
tive polarization. Given a source emitting a beam of Dirac 
particles whose polarization oscillates with a certain fre- 
quency f2, the different time of flight for particles with 
different spins gives a correction to the polarization mea- 
sured by a distant observer. For a beam moving in a Kerr 
space-time along a trajectory on the equatorial plane, we 
find 

Su) _ 3loL 2 Mci\ c 

~u 2^4 ' ^ ' 

M and a are the two parameters of the Kerr solution, 
the mass and the angular momentum per unit of mass, 
respectively, while r is the radius of the trajectory, and L 
is the distance traveled. For the Earth one gets 



*!L„ uL 2. 10 -4T_L 



(24) 



However, the realization of experiments on polarized 
beams of Dirac particles that are capable of testing such 
small deviations is rather remote (for photons similar de- 
viations have been confirmed experimentally, see [24] and 
references therein). 

VI- Concluding remarks. — In this work we have 
focused our attention on the semi-classical limit of the 
Dirac equation in a curved space-time. The main point 
is the way this limit has been taken. While the eikonal 
approximation is the standard tool to derive classical me- 
chanics from quantum mechanics, the hypothesis of local- 
ization is proper of our treatment, since other approaches 
do not consider it [5, 15]. But we want to stress that to 
make a comparison with the behavior of classical objects, 
a classical trajectory must come out from the quantum 
description. 

We obtain such a feature by a multi-pole expansion and, 
since we are describing an intrinsic spin, we have to retain 
only the monopole term when taking the classical limit. 
The best we can do is to localize wave functions at the 
Compton scale order, but we were able to define quanti- 
ties on their mean values by a spatial integration. Such a 
kind of integration is proper of the Papapetrou approach 
and this way a trajectory is well-defined. 



This formulation provided with a notion of spin tensor, 
whose quantum character is manifest, because it arises 
as a first order correction in the semi-classical approxi- 
mation, and whose physical interpretation is clear as the 
expectation value of the spin density. Furthermore, gen- 
eral covariance is manifest during all calculations, while no 
weak-field approximation takes place for the gravitational 
field. 

The dynamics reproduces the results by Mathisson and 
Papapetrou at the first order in the spin. Hence, there is 
no agreement with the papers [7] , where the semi-classical 
limit is performed by a Foldy-Wouthuyscn transformation 
on the Hamiltonian and the Pomeransky-Kriplovich for- 
mulation is obtained. These two different procedures al- 
low us to infer a semi-classical picture, but they provide 
different results for the classical dynamics. 

Therefore, it is possible to infer from a quantum descrip- 
tion both the two compelling formulations for spinning 
particles in GR. Only the experimental detection of devia- 
tions from geodesic motion would tell us what is the appro- 
priate semi-classical procedure for Dirac fields in curved 
space-time. 
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